Step-like discontinuities in Bose-Einstein condensates and Hawking radiation: 

dispersion effects 



O 



X 



Carlos MayoraQ 
Departamento de Fisica Teorica and IFIC, 

Universidad de Valencia- C SIC, 
C. Dr. Moliner 50, 46100 Burjassot, Spain 

Alessandro Fabbrfj] 
Departamento de Fisica Teorica and IFIC, 

Universidad de Valencia-CSIC, 
C. Dr. Moliner 50, 46100 Burjassot, Spain 
and 

APC (Astroparticules et Cosmologie), 
10 Rue A. Domon et L. Duquet, 
75205 Pans Cedex 13, France 



Massimiliano Rinaldlj] 
Department de Physique Theorique, 
Universite de Geneve, 
24 quai E. Ansermet CH-1211, Geneve 4, Switzerland 

(N I (Dated: May 30, 2011) 

In this paper we extend the hydrodynamic results of [l| and study, analytically, the propagation 
of Bogoliubov phonons on top of Bose-Einstein condensates with step-like discontinuities by taking 
into account dispersion effects. We focus on the Hawking signal in the density-density correlations 
5-H , in the formation of acoustic black hole-like configurations. 



I. INTRODUCTION 



> 

in 

The study of analog models of gravity in condensed matter systems ip, y| has motivated the investigation of quantum 
effects in gravity, in particular Hawking radiation from black holes [J] , in the presence of modified dispersion relations 
(see 0-0] 1^ Modified dispersion relations at high frequency have also been considered in many papers in cosmology 
(see e.g. Q), but also in the context of the Unruh effect [9|, which is closely related to the Hawking emission from 
a black hole. On a more formal level, issues related to quantum field renormalization in the presence of dispersion 
were investigated in [Tfjj . Among the many systems proposed to create black hole- like configurations, e.g. superfluid 
liquid Helium 11], atomic Bose-Einstein condensates (BECs) [Hj], surface waves in water tanks [ljj, degenerate Fermi 
gases (l4| . slow light in moving media (l5j . traveling refractive index interfaces in non linear optical media [l6j], BECs, 
characterized by superluminal dispersion relations, appear to be quite attractive from the experimental point of view 
[Tt| . In this context, recently an alternative measure of the Hawking effect was proposed in terms of non local density 
correlations [l8[ for the Hawking quanta and their partners situated on opposite sides with respect to the acoustic 
horizon. The calculations were performed using the gravitational analogy, which corresponds to the hydrodynamic 
approximation of the theory. This proposal was validated with numerical simulations within the microscopic theory 
(19 1, in dicating that the Hawking signal in the correlations is indeed robust. Subsequent investigations were performed 
in 2(J (where analytical approximations based on step-like discontinuities were considered) and [2l| using stationary 
configurations. 

In this paper we extend the hydrodynamical analysis in ll| and consider in particular the effects of the temporal 
formation of acoustic black hole-like configurations, as in [18[ and including dispersion effects. Our analytical 
analysis is based on step-like discontinuities in the speed of sound and thus extends the stationary results in [20j . We 
mention that step-like configurations in BECs were also considered in [2a - l26l |. 
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The plan of the paper is the following: in section II we briefly describe the model used and the basic equations, 
while in sections III and IV we analyze thoroughly the stationary case (spatial step-like discontinuities) and the 
homogeneous one (temporal step- like discontinuities). By combining the results of these two sections, in section V 
we discuss the main Hawking signal in correlations for the formation of acoustic black hole-like configurations and in 
section VI we end with comparisons with the hydrodynamical results in |18| . 



II. THE MODEL AND ITS BASIC EQUATIONS 



We start with the basic equations for a Bose gas in the dilute gas approximation described by a field operator $ 
[27]|-[29|]. The equal-time commutator is 



$(t,x),¥(t,x')] = 5 3 (x-x') 
and the time-dependent Schrodinger equation is given by 



ihdt^ = ( -— v 

2m 



(1) 



(2) 



where m is the mass of the atoms, V ex t the external potential and g the nonlinear atom-atom interaction constant. 
By considering the mean-field expansion 



* ~ * (1 



(3) 



with (f) a small perturbation. The macroscopic condensate is described by the classical wavefunction which satisfies 
the Gross-Pitaevski equation 



2m 



^ + Vext + gn * , 



(4) 



where n — \^q\ 2 is the number density, and the linear perturbation 4> satisfies the Bogoliubov-de Gennes equation 



«ft5 t = - — V 
\ 2m 



h 2 h 2 V* 



m 



V J 4> + mc 2 [ (4> + 4> f ) , 



(5) 



where c = is the speed of sound. 

V m ^ i — i 

To study analytically the solutions to ([3]), along the lines of [19|, we shall consider condensates of constant density 
n and velocity (for simplicity along one dimension, say x). Non-trivial configurations are still possible, provided one 
varies the coupling constant g (and therefore the speed of sound c) and the external potential but keep the sum 
gn + Vext constant. In this way, the plane- wave function "Jo = yfne lkaX ~ lw ° t , where v = is the condensate velocity, 
is a solution of Q everywhere. 

The non-hermitean operator </> is expanded as 



(G) 



where Oj and fit are the phonon's annihilation and creation operators. From ([5]) and its hermitean conjugate, we see 
that the modes 4>j(t,x) and (pj(t,x) satisfy the coupled differential equations 



i(dt + vd x ) + ^d 2 x -~ 



-i(d t + vd x ) 



^d 2 -- 
2 x tl 



TV j 



(7) 



where £ — h/(mc) is the so-called healing length of the condensate. The normalizations are fixed, via integration of 
the equal-time commutator obtained from ([T]), namely 



[4>{t,x),fi{t,x')\ = -S(x-x') , 



(8) 



3 



by 



/ 



dxfaQ, - $<p jf ] = . (9) 



We shall consider step-like discontinuities in the speed of sound c, which is the only non-trivial parameter in this 
formalism, and impose the appropriate boundary conditions for the modes that are solutions to Eqs. ([7]). A similar 
analysis was carried out in the hydrodynamic limit £ — ► in the work [l|, by using the more appropriate density 
phase representation 

? n 1 1 

*=Tn +l lf (10) 



III. STEP-LIKE SPATIAL DISCONTINUITIES (STATIONARY CASE) 

In this section we study dispersion effects for the case of spatial step-like discontinuities. We treat subsonic configura- 
tions in subsection IIIIA1 thus extending the hydrodynamic analysis of [l|, and subsonic-supersonic ones in subsection 
IIIIBI This case is particularly interesting in view of our application to study the main Hawking signal in correlations 
from acoustic black holes, along the lines of [l9j . 



A. Subsonic configurations 



We consider a surface (that we put for simplicity at x — 0) separating two semi-infinite homogeneous condensates 
with different sound speeds: c(x) — ci0(—x) + c r 9(x). The velocity of the condensate is taken to be negative (v < 0), 
so that the flow is from right to left. We assume that the condensate is everywhere subsonic, that is \v\ < c r m, and 
that v, ci and c r are time-independent. 

To explicitly write down the decomposition of the field operator 0, we first need to study the propagation of the 
modes and construct the "in" and "out" basis. To understand the details of modes propagation, we need to solve 
the equations ([7]) in the left and right homogeneous regions, and then impose the appropriate boundary conditions. 
These simply are the requirement that <fi and ip, along with their first spatial derivatives, are continuous across the 
discontinuity at x = 0. 

We denote the modes solutions in each homogeneous region and corresponding to the fields <j> and tp as ]j e - lwt + ikx 
and E e ~ lwt + lkx respectively. The boundary conditions at the discontinuity, as we will see explicitly later, require us 
to work at fixed uj. Therefore we write the modes as 



so that the equations ([7]) simplify to 



(11) 



(w — vk) 
(w — vk) 



£ck 2 
^ck 2 



while the normalization condition (J9j) (J = uS] gives 

\D(u)\ 2 -\E(u)\ 2 



D(u) 



E(w) = -£>(w) 



I 



2ntvn 



dk 



dw 



The combination of the two Eqs. (|12[) gives the non linear dispersion relation 

/ f 2 k A " 

I i \2 2/12,?™ 

[w — vk) 



(12) 



(13) 



(14) 



plotted in Fig. [TJ At low momenta (k <§C |) we recover the linear relativistic dispersion, while at large momenta 
(k ^> |) the nonlinear superluminal term becomes dominant. Moreover, inserting the relation between D and E from 



4 



3 o.o 
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FIG. 1: Dispersion relation for subsonic configurations. The solid (dashed) line corresponds to the positive (negative) norm 
branch: u — vk = +{—)cd k 2 + 



12]) into (|T3f we find the mode normalizations 

£>(«) = ■ 

E(u) = - 



uj — vk + 



A-Khnc^k 2 (aj-vk)(^) 1 



w-vk-&- 



Anhnc^k 2 



(15) 



where k = k(w) are the roots of the quartic equation (fT4"]) at fixed w. Eq. (fT4")l admits, in the subsonic case, two real 
and two complex solutions. Regarding the real solutions, we will call k v and k u the ones corresponding to negative and 
positive group velocity v g — ^ respectively. They admit a perturbative expansion in the dimensionless parameter 



namely 



k v 
k it 



UJ 



v — c 
V + c 



c 3 z 2 



(v — c) 3 

(?Z 2 

(v + c) 1 



+ 0(z 4 ) 
+ 0{z A ) 



(16) 



The other two solutions are complex conjugates. We call kd(k g ) the roots with positive (negative) imaginary part, 
which represent a decaying(growing) mode on the positive x > axis and a growing(decaying) mode in the negative 
(x < 0) one. Such roots are non-perturbative in £ as they diverge in the hydrodynamic limit £ = 0, when Eq. (|14p 
becomes quadratic. However, they admit the expansions 



k 



d{g) 



(c 2 + v 2 )c*z 2 
4(c 2 - v 2 ) 3 



+ 0(z 4 ) 



+ (-) 



2i v / c 2 " 



(c 2 +2t> 2 )c 4 z 2 

Of„2 _ ,,2\3 



+ 0(z 4 ) 



(17) 



In what follows we do not need to specify the normalization coefficients for these modes, that we call generically 
^ v) h and ^ iv) h for the decaying and growing modes respectively, of the fields cf> and (p. 
In summary, the most general decompositions of <f) and (p in the left and right regions are given by 



l(r) 

UJ 

(r) 



(18) 
(19) 
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The coefficients A£ d G are the amplitudes of the modes, not to be confused with the normalizations coefficients. 
Indeed, the latter are determined uniquely by the commutation relations and the equations of motion, while the 
amplitudes depend on the particular choice of basis, as shown below. The matching conditions at x — to be 
imposed on Eqs. ([7]) are 



[</>] = o, 



0, [ip] = 0, [ip'] = 0, 



(20) 



where [ ] indicates the variation across the jump. It is clear that these conditions require w to be the same in the I 
and r regions. Eqs. (|20j) can be written in matrix form 



where 



IK 

K 

\A l d J 



= W r 




(21) 



Wi = 



( D ! v l 








ik v D v 








K 


EL 




< 




<e 1 u 







(22) 



and 



W r = 



I 


D r 




a 4> 


% \ 




ik r v D r v 












E r u 


ik r d cT v 




\ 


iky El 


iW 



(23) 



Multiplying both sides by Wr 1 we have 



A l u 
A'c 
\A l d J 



= M, 



( K 

K 

A r 
\ d 
\Al 



(24) 



The 4x4 matrix M scatt = W l ~ 1 W r encodes all non-trivial scattering effects due to the matching conditions (|2"0"|) . The 
form of M aca tt is much more involved than that found in the hydrodynamic limit in |l| . 

To construct M scat t we have used the general decompositions (|T5)) and (fr9|) . Not all modes, however, are physically 
meaningful. The validity of the mean-field approximation Q implies that only spatially bounded modes have to be 
taken into account. This means that the amplitudes of the growing modes (that diverge exponentially in the I or r 
regions) must be set to zero. There are no constraints, instead, for the amplitudes of the decaying modes. Indeed, as 
we will see explicitly in the construction of the "in" and "out" modes basis that follows, by taking into account the 
(I and r) decaying modes we have each time four amplitudes which are uniquely determined by our four matching 
equations. The physical meaning of the decaying modes is to "dress" the "in" and "out" modes basis, and this affects 
the calculation of local observables (this discussion follows that of [30j). 

We now proceed to construct the "in" and "out" modes basis for the case v = in a perturbative expansion up to 
0(z 2 ). This case can also be treated exactly, as shown in the Appendix A. The perturbative construction of the "in" 
modes for the more complicated case v ^ is given in Appendix B. To appreciate similarities and differences with 
respect to the hydrodynamical case treated in [l|, let us construct perturbatively the "in" and "out" modes basis, 
displayed schematically in Fig. [5] We consider the modes of the field <j). An identical analysis is valid for ip, up to the 
replacement of the D —> E. 

v,in 
1 u>,<j> 



The 

Die 



m 

-iujt+iklx 



w-mode u', is defined by an initial unit-amplitude left-moving w-mode coming from the right (= u^. 



), which is partially transmitted into a u-mode in the left region (u^ = D l v e 



) with amplitude A l v 



and partially reflected into a right-moving u-mode {u^'l — Die ZUJt + l K x "j with amplitude A r u . The construction is not 
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FIG. 2: "in" and "out" modes for spatial step-like discontinuities between homogeneous subsonic regions. We display the 
propagating modes (straight lines) and the decaying modes (curved lines), along with their amplitudes. 



finished yet, as we need to include as well the decaying modes in the left and right regions {u^f^ — D r P^e iut+tk d( g ) x ) 
along with their amplitudes A r d and A l d . In this way we have a total of four amplitudes which are uniquely determined 
by solving the following system of four equations 








M.. 



( 1 

K 



(25) 



V o 

By treating M scaU perturbatively in the parameter zi = we find, up to 0(z 2 ), the following solutions 

iy/ci (ci - c r ) 2 zi ci (ci - c r ) 2 (cf + cf) zf _ 



A ( . 
A 1 , 



_ T 

Cl+Cr C 3/ 2 ( C;+Cj ,) 2c 3 r (ci+C r ) 2 

Q - c r ici (q - c r ) 2 zi ci (q - c r ) (2cf - 3cfc r + 2qc 2 -I- cf) zf 

Ci + C r C 2 (q + C r ) 

(Q - C r ) ^fzi (Cl - C r ) zf 



d^V^ ( C l + °r) Id^cJ 2 (Cl + C r ) 
Cl (-Cl + C r ) yfzi cf (Cl - C r ) zf 



4cf (Cl + C r ) 
[c,. 2 + i (c 2 + C 7 2 - 



R 



c r ci)] = D L , 



dLcJ 2 (Cl + C r ) 2d\cJ 2 (Cl + C r ) 



[q +i(ci - 2c r )} = D R > 



(26) 

(27) 
(28) 

(29) 



In the limit z\ — > 0, we recover the results of As we can see, the amplitudes of the asymptotic modes A l v and 
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A T U develop an imaginary O(zi) contribution plus a real O(zf) one. These combine in such a way that the unitarity 
relation \A l v \ 2 + \A r u \ 2 = |i?| 2 + \T\ 2 = 1 is satisfied non trivially at O(zf), as 



a i2 



4c ; c r | u 2 (ci-c r ) 2 (c 2 + c 2 )tf 
2c ; c r (c; + c r ) 2 
^ 2 ( Q -c r ) 2 (c 2 + c 2 )g 2 

2c ;C 3 (q + Cr) 2 



(C; + C r ) 2 



Q + C r 



(30) 
(31) 



Finally, note that, although the amplitudes of the decaying modes do not enter in the unitarity relation, they are part 
of the full mode and give contributions, for instance, in the computation of density-density correlations. 

. Mode t#j 



m (/-iixxip // j, is composed by an initial unit-amplitude right-moving u-mode (u^ = D l u e vat + % K x ) coming 



The "in" u-mode 

from the left, along with the transmitted u-mode (u% r ) with amplitude A r u and the reflected u-mode (u^J 1 ^) with 

amplitude A l v . Here too we have decaying modes, with amplitudes A d ,A l d . All these amplitudes are obtained by 
solving 



1 





AT, 




and, up to O(zf), we have 



(32) 



A 1 
A r 

A 1 

a:, 



c r - q _ i{ci-c r ) 2 zi + (Cj - Cr) (cf + 2c 2 c r - 3c ; c 2 + 2cj?) zf 
Q + c r c r (ci + c r ) 



4 C 3 (q + Cr) 



R' 



2y/clc; . 7q (ci - c r ) 2 zi sjcx (q - Cr) 2 (cf - 4c;c r + c 2 ) zf = 



Q + c r 



3/2, . \ 
C r ' (C; + C r ) 



+ ■ 



(q - C r ) 



(C; - C r ) 

^x/Q (cj + Cr) ' 2d l <t> y/cic r (ci + C r ) 
^fc~l(-Ci+ C r ) ^fz[ y/ci(ci-Cr) 



8cl /2 (Cl + C r ) 

{-c r + i (2q - c-)] zf = D L , 

[q 2 + i (q 2 + C r 2 - QC r )] zf = D R . 



(33) 
(34) 
(35) 
(36) 



d^Cr (Cl + C r ) 2e^q! (a + Cr) 

The unitarity condition for the asymptotic modes \A l v \ 2 + |^| 2 = |i?'| 2 + \T'\ 2 = 1 is again non-trivially satisfied, as 

' C r - Ci \ 2 Cl (Cl - C r ) 2 (c 2 + C 2 ) zf 



\Ar = 



AH 2 = 



Cl + C r J 2cl (ci + C r ) 2 

AciCr Cl (C X - C r ) 2 (c 2 + C 2 ) Z 2 



(Cl + Crf 



24 (Cl +C r ) : 



(37) 
(38) 



. Mode u v 'T 



The "out" u-mode u v J is made of a linear combination of initial right-moving (u"'^) and left-moving (u v J r ^) com 



u,4>> 



ponents, with amplitudes A u and A T V , producing a final left-moving w-component (u v J ^) of unit amplitude. The 
amplitudes, together with those of the associated decaying modes, are given by solving 



( \\ 

A 1 

u 



\ Al J 



= Af, 




(39) 
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and, at 0(z 2 ), one has 

A 1 - 
u — 

A r - - 
A 1 - 

a:, = 



c r - Q i (ci - c r ) 2 zi (c; - c r ) (cf + 2c 2 c r - 3qc 2 + 2c^) z 2 = rU 

Cl + C r C r (Q + C r ) 4c^ (c| + C r ) 

2^qc7 i^Q (q - C r ) 2 Zi (Cj - Cr) 2 (c 2 - 4c ( c r + c 2 ) z 2 _ , 



Cl + c r 



Cr /2 (C; + C r ) 



%cj 2 (c; + C r ) 



+ 



[c 2 - i (q 2 + a 2 - c r ci)] = D* R . 



d^Cr (C t + C r ) ' 2d^ C 3 (a + Cr) 

One can easily check that the unitarity relation is satisfied, as \A l u \ 2 + \A^\ 2 = \R\ 2 + \T\ 2 = 1. 



(40) 
(41) 
(42) 
(43) 



Mode u 



u,out 



We finally consider the mode u 



u,out 



This is defined by initial right-moving and left-moving components, with 



amplitudes A l u and A r vl resulting now in a final right-moving u component (u"'^) of unit amplitude. The system of 
equations to be solved (taking into account the decaying modes) is 



A u 


\ A d) 



M 



scatt 



( K 

\ 



and its solutions, up to 0(z 2 ), are 



(44) 



A 1 - 

A v — 

A 



2y/Qc7 | iy/ci(ci - C r ) 2 Zl y^Q (cj - C r ) 2 (cf - 4c;C r + cf) zf _ 



+ 



C l + C r cj 2 (CJ + Cr) 8cl /2 ( Q + C r ) 

ci-Cr ici(c t - c r ) 2 zi Q (cj - c r ) (2cf - 3c 2 c r + 2qc 2 + cf) zf 



d Ji 



Cl + C r C 2 (Cl + C r ) 

Cl (ci - C r ) Zl , (Cl - C r ) zf 



R* 



ly/cT (Cl + C r ) 2d l 4> C r /2 ( Cl + C r ) 
Ci (-Cl + C r ) Zi cf (a - C r ) zf 



d^Cr (C| + C r ) 2d r ,cV 2 (Cl + C r ) 



Acj. (q + C r ) 
\-c 2 + i (c 2 + c 2 - cicf)] = D L >, , 

[a + i (2c r - ci)] = D rU , 



(45) 

(46) 
(47) 

(48) 



with unitarity condition \A l u \ 2 + \A r v \ 2 = \R\ 2 + \T\ 2 = 1 satisfied up to 0{zf). 

Having constructed explicitly the complete "in" and "out" modes basis, we can now write the two alternative 
decompositions for the field operator <j> 



r-OO 

4> = / diO hvJn(out) u v,in(out) ^ ^ + & u,in(out) u u,in(out) ^ ^ + 
JO L 

_|_ ^v, in (out) An (out) *^ ^ ^_ ^u,in(out)^ u u,in(out)* ^ ^ 



(49) 



The relations between the "in" and "out" modes are 



u. 



u. 



rj-\ V.OUt . T-> U.OUt 



Tl u,out 

LJ.<p ' 



71/ V.OUt 

R u 



and are valid for all components of the modes basis, decaying modes included. This allows us to find 

~v,out rp~v.in | 7V ~uAn 

~u,out r>^.u,in | rp/^u.in 

a = Ra + 1 a,' . 



(50) 



(51) 
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Density-density correlations 



The basic quantity that we want to study in detail later is the one-time, normalized, symmetric, two-point function 
of the density fluctuation 



G (2) (t;x,x') = limHIri^^nV^Ollin) 
2n z t-tt' 



(52) 



where {, } denotes the anticommutator, and the the operator n 1 = n((f> + (f)') (see eq. (jTUJ) ) can be expanded in the 
two equivalent "in" and "out" representations, 



h 1 = n 



did 



~v,in(out)( v,in(oui) v,in(out)\ , ~u,in(out) j u,in{out) uSn(out)\ , v, _ 



Thus, the general two-point function in (|5"2"j) explicitly reads 



{in\{n 1 (t,x),n 1 (t r ,x')}\in)\ = it 



did 



where 



n, V ' m _|_ 1 , V ' m 



u,m . u,in 



= e - *"* \(D r v + E r v )e lk ^ )x + R{D r u + E r u )e lk ^ u)x + T{D l v + E l v )e lk ^ )x + 

+ {Di,d\ + Di,d l v y k >^ + (DUl + D%dr v )f&w*_ , 

= e - *"* \{D l u + E l u )e< {u)x + R'(D l v + E l v )e tk ^ u]x + T\D r u + E r u )e ik ^ x + 

+ (D*^ + Dlfye^W* + (D% + D*d r v y k ^ x . 



(53) 



(54) 



(55) 



Let us consider, for instance, one point located in the left (a; < 0) region and one in the right (x 1 > 0) one. Substituting 
the expressions above into (154[) . we see that there are u — u and v — v contributions, while the u — v term, being 
proportional to R*T + R'T'*, vanishes. Finally, the contribution coming from the decaying modes is subdominant. 
Therefore, the integral (|54j) is well approximated by the hydrodynamic approximation, obtained for small uj, namely 



G {2) {t;x,x') ~ - 



2irmn(c r + q) 



(V - Cl )(v - Cr) + ^) (V + Ct)(v + c r ) (-^ + ^-J 



(56) 



B. Subsonic-supersonic configuration 



Unlike the spatial step-like discontinuities studied in [lj in the hydrodynamical limit, dispersion effects allow us to 
study also configurations with supersonic regions. Since we are interested to model black hole-like systems, we shall 
consider the case where there are one subsonic and one supersonic region separated by a sharp jump in the speed of 
sound. Therefore we write c{x) = ci6{—x) + c r 9(x), where now q < |u| and c r > \v\. The modes in the subsonic 
region (x > 0) are the same as in the previous subsection. In the supersonic (x < 0) part the dispersion relation (fT4")) 
changes and it is represented in Fig. [3] We see that, for u> less than a certain value that we call u) max , there are now 
four real solutions, corresponding to four propagating modes. Two of them are present also in the hydrodynamical 



approximation, and, when expressed through the variable zi 



6" 



, they read (we omit the subscript I) 



h v — 



v — c 

UJ 

v + c 



c 3 z 2 
(v + c) ; 



+ 0(zf) 



+ 0(zf) 



(57) 
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FIG. 3: Dispersion relation in the supersonic case. Positive (negative) norm modes belong to the solid (dashed) line. 



and, unlike in the subsonic case, they both move to the left, as 4r \k u(v) < 0. The value k v belongs to the positive 
norm branch while k u to the negative norm one, as shown in Fig. [3J The other two values of fc, called and k4, exist 
because of dispersion, and are not perturbative in £. In fact 



«3(4) 



1 



(c 2 + v 2 )c A z 2 
4(c 2 — v 2 ) 3 



+ 0(z 4 ) 



+ (-) 



2Vv 2 - c 2 



1 



(c 2 + 2v 2 )c i z 2 
8(c 2 - v 2 ) 3 



+ 0(z i ) 



Comparing with the expressions (|17[) , we see that k% and k^ are the analytic continuation for supersonic flows of the 
decaying and growing modes seen in the subsonic regime. These two modes (which belong respectively to the positive 
and negative norm branches of Fig [3]) both move to the right as ^U 3(4) > 0. This means that they are supersonic 



and able to propagate upstream, against the direction of the flow. The value of w r , 
explicitly by imposing 4r:|fe=fc max = 0, where 



U)(k 

max 

) can be calculated 



v 2 v v 2 

-2 + — \/8+ — 



2c 2 



2c 



1/2 



(58) 



One can easily check that u> max and k max are well inside the non-perturbative region (~ l/£). When lu > uj maxi 
instead, we find again two real propagating modes (k real) and two complex conjugate ones, corresponding to decaying 
and growing modes, just like in the subsonic case. Thus, for u > uj max the analysis is the same as in the subsonic 
case, and so we omit it. 

Let us now write the general solutions for cj) and tp in the left (I) and in the right (r) regions for lj < uj max . In the 
Z-region we have 



D[A l v e ik * 



E l v A l v e ik * 



D l A 1 e" 



E l u A l u e ik « 



+ D^A 1 ^ 1 



ElAle ik 3 x + E l d A 



DiAie ik * 



l Al 



while in the r-region we find 



D r v A r v e lk > 



E„A„e « 



D r u A r u e ik > 



ElA T J K 



- d 4 A r d e ik * x 
d v A d e lk > - 



G^A r g e ik > 



G v A r g e l 



The D and E normalization coefficients of the propagating modes (four in the supersonic region and two in the 
subsonic region) are given by Eqs. (fT"5|) . As before, the matching conditions (|20|) can be written in the matrix form 



Wi 



( A \\ 
A 1 

4 



= W r 




(59) 
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where W r is the same as Eq. (|23|) . while Wi is given by 



I 


D l 


D l u 








ik v D v 


ik l u D l u 


ik 3 D 3 






■ f% 


EL 


■f\ 


E{ 


V 


ik v E v 


<El 


ik a E 3 


ik\E\ J 



(60) 



Multiplying both sides by W, 1 we find 



A 1 



= M„. 



( K 

A r 



(61) 



where M scaU = W l ~ 1 W r encodes the scattering effects due to the matching conditions ([20]) . As in the previous 
subsection, we shall proceed to the construction of the "in" and "out" mode basis for this configuration. With these, 
we will construct the decompositions of the field 4> along with the the density-density correlations. 



Construction of the "in" and "out" basis 



We shall now construct the "in" and "out" basis, which are now composed of three modes each, as shown in Fig. 2] 
Below, we find the leading-order amplitudes of the various amplitudes. In Appendix C, we display the next-to-leading 
order terms for it^'T and MjT* in order to show that unitarity relations are non-trivially recovered. 




FIG. 4: 'in' and 'out' basis in the subsonic-supersonic configuration. 
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Mode <7 



The mode u^™ is defined by an initial left-moving unit-amplitude component (u^) coming from the subsonic region 



v,r ■ 



on the right, which generates a reflected right-moving mode (u™'^) with amplitude A r u , together with the associated 

decaying mode with amplitude A r d . In addition, now there are two transmitted modes, one with positive norm (u^ 1 ^) 

and the other with negative norm (u"'^ ), with amplitudes A l v and A l u respectively. These can be computed by solving 
the system of equations 



A 1 


V / 

The leading order 0(1) solution in a zi expansion is 
fcTv-Ci 

A r 
A 1 

a: 



m, 



scatt 



t 1 

K 

A r d 

V o 



(62) 



Q V — C r 

v + c r 

v — c r 

ACr V + Q 
Q C r — V 



ciJz~i^c r {v 2 - cf) 



v+^Jv 2 -cf 



^ («-Q)( C 2-i; 2 )3/2( Cr+Q 

The amplitudes of the propagating modes satisfy the unitarity condition \A l v 



+ i [v^jv 2 -cf + v 2 - c\ 

Kl 2 = i- 



(63) 



\a: 



• Mode u 3 '™ 

UJ,(p 



The mode 



3, in 
l uj,<t> 



corresponds to a unit amplitude, supersonic positive norm right-moving plane wave from the left 



(uf,'^ = D\e «"*+ lfc 3(") a; ) j which is reflected into a positive norm (u^J 1 ^) and a negative norm (u"'^) component with 
amplitudes A% and A^ moving to the left. In addition, there is a transmitted right moving mode in the subsonic 
region (w"'^,) with amplitude A 1 ^ and the decaying mode with amplitude A r d . By solving the system 



(64) 



(Ai\ 

A L U 
1 




( \ 


— M scat t 


AR 

A R d 


V J 




V / 



we find, at leading order in z\, 



{v 2 -c 2 fl A {v + c r ) 



AR - 



A,, = 



Af = 



c 3/2 \Z2z~i(ci + c r )yjcl - i 

^2F r {v 2 -cfY^jv + Cr) 
Cl^fziic 2 - C 2 )y/c 2 - V 2 

(v 2 -cf) 3 / 4 (v + c r ) 



v 2 + i 



f 2 



c^ 2 ^/2z~i(ci - 
(v 2 - c 2 ) 1 ^ 
2d^{v 2 -c 2 ) 



\Jc 2 -v 2 +i 



l f 2 



(v - iy/c 2 . - v 2 ). 



(65) 



Note that the amplitudes of the propagating modes diverge in the zi — > limit, and that, at leading order in zi, one 
has \A^\ 2 + \A^\ 2 - \A^\ 2 = 0. In order to check the unitarity condition \A l v \ 2 + \A r u \ 2 - \A l u \ 2 = 1, we need the 
next-to-leading order expansion, which is displayed in Appendix C. 
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Mode u 



4,in* 
w,4> 



The mode 



4,m* 
l u,4> 



(where * means that this is a negative norm mode) consists of an initial unit amplitude supersonic 
right-moving component from the left (it*'^ = D l 4 e~ uot+ ' lk *^ x ) , generating a reflected positive left-moving norm 
mode {u v J 1 ^) and negative norm left-moving mode (u™'^) with amplitudes A l v and A l u respectively. Moreover, in 
the subsonic region one has a transmitted right-moving wave (it"'^) with amplitude A r u , and a decaying mode with 
amplitude A r d . By solving 



we find, at leading order, 



a 1 : 



M 



scatt 



A 1 ' 


V i / 



{v 2 -c 2 f/\v + c r ) 
c^ /2 ^/2zi(ci + c r )^/c 2 r - v 2 
V2E;.(v 2 ~c 2 f' A {v + c r ) 

Cly/zi(Cr ~ C 2 )V c r ~ y2 

{v 2 -c 2 f/ A {v + c r ) 

3/2 



V o / 



i^Jv 2 



V4 



i\J~v 2 



(v 2 - c 2 fl A {v 2 -c 2 + v^v 2 - cf) 
2d 4> {c 2 . - v 2 )(c 2 -v 2 + v^F^c\) 



(v — i\J c 2 — v 2 ) 



(66) 



(67) 



As for 



3, in 



the amplitudes of the propagating modes diverge when zi 



\a: 



and at this level of approximation they 
\A l u \ 2 — —1 is checked in the Appendix C 



satisfy \A l v \ 2 + \A r u \ 2 - \A l u \ 2 = 0. The unitarity condition \A l v \ 2 
by considering the next to leading order terms. 

The construction of the "out" modes proceeds similarly. These are u^'^J 4 ', u^^ ut (of positive norm) and '^ ut * (of 
negative norm), which are composed by appropriate combinations of initial right- moving and left-moving components 

u.l* 



(plus the associated decaying mode). These generate respectively, unit amplitudes u^J^, u™'^. 
detail, we have the following cases. 



and 



More in 



Mode u 



u).<t> 



In this case, one needs to solve the system 



( 1 \ 



4 



M, 



scatt 




(68) 



which yields, at leading order, 

A 1 -- 

A 1 - 
A r — 

a:, = 



(v 2 - eff^ ^cj^ 

v)(Cr + Cl 



2zic l /2 (c r 



(v 2 



cf) = A[ 



2zic l /2 (c r - v)(c r + ci) 
~c7ci-v _ i 

— A v > 

Cl c r — V 



i\fcl 



v 2 ) = iAt 



(v 2 - cf)(c 2 . - V 2 + Vyjv 2 - IJ) 

V2d<f,(cr - v)(ci + c r )y/ci(v 2 - c 2 )ci^/Ti 



(69) 
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• Mode u 



In this case the system to solve is 



and the solutions are 



A 1 - 

A 3 — 



4' 



A r v — 

a: = 



( ° \ 

K 
o 

\ A U 



(K 
l 

V 



/2^(« 2 - cf ) 3 / 4 (« + cv) 



2b~ r {v 2 - C 2) 3 / 4 (t; + c r ) 



— l\ v — c 



2\ - 



— A,. 



yfz[ci{c 2 
v + Cr 



(-V4 



i^v 2 -c 2 ) 



A' 



, V2c r (v 2 - c 2 )(c 2 r -v 2 + v^v 2 - c 2 ) 
d^Cr - v) 3 / 2 (c 2 - c^) v /c r (w + c r )ciy/z~i 



(70) 



(71) 



• Mode u 



ul, out* 

UJ,(p 



In this final case, the system is 



( \ 
i 

A{ 



( K 
o 

V 



and the solutions read 



A< 



A = 



2zic,' (v - c r )(c r - q) 



-A., 



f2zic{ {v-c r ){ 
{v 2 -c 2 fl^ 2 



Al 



A r — 



2zic l /2 (v - c r )(c r - q) 



C r Cl + V 
Cl V - C r 
2 „,2 



(-V4 



-A 



(cf — V 2 )(V 2 — C 2 — Vy/v 2 — C 2 ) 



VZd^v - c r )(ci - c r )^ci(v 2 - c 2 )c iy /Ti 
With these results, we are able to write down the relations between the "in" and "out" modes 



n v,out , r ur,out 



a\ ul.OUt* 



3,m 



4,m* 



A L v. out , a R ur.out i a L ul.out* 



u,4> 

A I' V,OUt 



at' ur,out 
A u 



V ul,out* 



(72) 



(73) 



(74) 



We note that, unlike the subsonic case (|50[) . we now have combinations of both positive and negative norm modes. 
Because of this, the two decompositions (we restrict our analysis to the case u> < uj max because it is the relevant one 
for our subsequent discussion) are given by 



duj 
dui 



~v,in v,in -3,m 3,in -4,m 4,in -»,inf »,in* , ~3,mt,,3irw , ~4,inf 4,int 
^.v, out v, out . <-ur,out ur,out , ~ul,out ul,out . ~v,outj v,out* , ~ur,outf ur,out* , ~ul,out\ ul,out* 



(75) 



(76) 
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and they are inequivalent. This can be easily seen by using (|74[) to find the relation between the two families of a and 
a' operators 



a 



^ur.out 
ul,out\ 



= A 1 a v ' m - 



A^al m + Alai m \ 

■ A^r + A'!aff nt . 



(77) 



The fact that the r.h.s. of these relations contain both creation and annihilation operators makes it clear that the 
two decompositions do not share the same vacuum state (|in) ^ |out)). 

Density-density correlations 

To compute the normalized density-density correlation analogous to Eq. (f52|) . we first expand the operator n 1 in the 
"out" decomposition 



h (t, x) = n 



zv,out( v,out , v,out\ i ~ur,out( ur,out ur,out\ , ~ul,out( ul,out ul,out\ , . 



(78) 



and we use the relation between the "in" and "out" operators (|77p. This gives the following two-point function in the 
| in) state 



cfcj ■ 



where, explicitly, 



v,out . v,out 

U OJ,<j> + U U,Lp 



ur,out . ur,out 



uLout* 



ul.out* 



A£(D£ + E^e^O")* + A l 4 (D l 4 + E l 4 )e ik '^ x + A r d (d^ + dQe ik ^ u > 

e~ wt [(£>; + ^)e lfe - (w)a; + A r v (D r v + E r v )e lk ^ x + 

A\{D\ + ^Je'*^) 1 + AiC-D^ + E{)e ik ^ x + A r d (d r ^ + d£)e**S(<")« 

A l 3 (D l 3 + E l s )e ik ^ x + A l 4 (D l 4 + E l 4 )e tk ^ )x + A r d {d r ^ + d r v )e ik ^ u > 



(79) 



(80) 



(81) 



(82) 



The coefficients A l v , A r v ,A l 3 , A\ and A r d are given, respectively, in ([7T|) and (|75|) . The analysis of the main 

correlation signals has already been performed in (20j |. We are interested in the correlation between w" ,r and w^,''* 
because this represents the main signal due to the Hawking effect (correlation between the Hawking quanta and their 
partners). We take x (x 1 ) in the left (right) region and evaluate the following integral 



{in\{n\t,x),h\t',x')}\in)\{uY o u u J*) = n 



'a 

+ (<AT + A^ u AT){u u J i : + u^*)(t, x)(u u J r ; + uZM, x') + c.c. 



(83) 
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The values of the above amplitudes are given in (l63l) . (|65|) and (j67|) . We also take into account that 

(1 =b. 4'* 4 r -I- A L * A R — 4''* 4 r ' — fl 



r ur,out ul,out~\ 



(84) 



where we have used the relation between the "in" and "out" operators given in (1771) . The term A l *A r u is subleading 
with respect to the other two terms, which go as 0(1/oj), given that the main contribution to the integral above is 
valid for small uj. Note also that the products A^*A R (and A l u *A r u ) are real at leading order. Therefore we have 



(in\{h\t,x),n\t\x')}\m)\( U y <^ u% 1 *) ~ An 2 / dco 

Jo 

and, at equal times, the normalized two-point function is 

1 



u *A r u Re 



G^(t;x,x')(uy ^u^ 1 *) 



(v 2 - c 2 fl 2 



[U , + u 



Sill 



( x x \ 

c V v-\-c. r V-\-Ci / 



Ann ci(v + ci)(v - c r )(c r - ci) 



(86) 



v+c r 



v+ci 



This result, which coincides with the one given in [20j . gives an estimate of the Hawking signal in correlations only for 
stationary configurations. Our aim is to perform a similar construction, but for acoustic black hole-like configurations 
which are formed at some time along the lines of the numerical analysis presented in [19J. 



IV. STEP-LIKE DISCONTINUITIES IN t (HOMOGENOUS CASE) 

In this section, we study correlation functions in the case of temporally formed step-like discontinuities between 
homogeneous condensates. In subsection IIV Al we consider condensates which remain subsonic at all times. In 
subsection II V Bl we turn to the more relevant case when the final condensate is supersonic. 



A. Subsonic configurations 



We consider a step-like discontinuity in t (say, at t = 0), separating two infinite homogeneous condensates: c(t) = 
Ci n 0(—t) + c ou tO{t). In this section we consider |^| < c in ^ out j so that the condensate is subsonic at all times. The aim 
is to determine the mode propagation at all times, and to define the "in" and "out" mode basis. The appropriate 
decompositions of our field (/) will be given afterwards. 

The general solutions in the "in" (t < 0) and "out" (t > 0) regions describing the fields </> and ip are of the form 
£) e -iwt+%kx an£ j j^ e -iwt+ikx _ r-Q^ bo^jaj-y conditions at t = require us to work at fixed k. Therefore we write 



4> k = D{k)e~ lw{k)t+lkx , tp k = E{k)e~ tw{k)t+ikx 



(87) 



for which Eqs. (J7J become 



while the normalization condition 



(oj — vk) + 
{uj — vk) + 
yields 



c£fc 2 
cik 2 



D(k) = --E(k) 



E(k) 



-Mk) , 



\D{k)\ 2 -\E{k){ 



1 



2-Khn 



(89) 



The combination of Eqs. 
normalization coefficients 



gives rise to the non- linear dispersion relation (|14[) represented in Fig 1, and to the 

u — vk 



D(k) 



E(k) 



1 



^ Amhnc^k 2 \ {uj — vk)\ 
y/ Anhnc^k 2 \ (lu — vk)\ 



(90) 
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Here, ui = ui(k) corresponds to the two real solutions to Eq. (|14l) . which is quadratic in to at fixed k. These read 



c 2 k 4 £ 2 

uj+(k) = vk + \lc 2 k 2 -\ --2- 



= vk - \\c 2 k 2 + , (91) 

where uj+ (k) corresponds to the positive norm branch, and w_ (k) to the negative norm one. Note that there are no 
normalizable mode solutions with complex k, because in the infinite homogeneous "in" and "out" regions they would 
correspond to modes which decay on one side but grow without bound on the other. Therefore, at fixed fc, the general 
decompositions of <f> and tp in the "out" and "in" regions are 



out(in) _ ikx 
Vk — e 



Dt ut{in) (k)A out(in) e-^T Hin H^ + D^WB^e-*-'™™*] , (92) 
Kut^mout^e- 1 ^^ + E~ u t ^ n) {k)B out{in) e-^ in) ^] . (93) 



For k > (< 0) we have a positive norm right-moving (left-moving) mode (lu — uj + (k)) and a negative norm left- 
moving (right-moving) one (uj — w_(fc)). According to (J7J, the matching conditions at t — are 



which can be written in matrix form 



where 



V out(in) 

(k) 1 

out(iri) 

(k) 

Multiplying both sides by W~ u \ we find 



[4>] = o, M = , (94) 



w ~* ' bZ )= w 4 t ) > ^ 



/ D + ... Jk) D~ ... Jk) \ 

rir I out (in) \ ' outlin) v / 1 fc\£i\ 

W out(m) = + ^ > . p _ i . . (96) 



ou t 
out 



(97) 



Explicitly, the Bogoliubov matrix M\> og = W J t Wi n reads 



71 t ( ~t~ ^out Qin ^out \ /no\ 

Mb ° 9 = 2^i m n out { n m - n out n m + n out ) - ( 98 ) 

where we define J7 otlt ( m ) = \ [UJ out ( m ) _ p or v — o we recover the formulas given in [26j]. 
Connecting the "in" and "out" basis 

The "in" and "out" modes basis are easily identified in terms of positive-frequency "in" and "out" modes (%" ( j ) OM *' ) = 

D tn(out) < y k ) e ~ lU+{k)t+lkx '> for u Zlp Ut) the analysis is identical up to the replacement of D+ n(out) (k) by E+ n(out) (k)) 
which are, respectively left-moving (fc < 0) and right moving (k > 0). To connect them, as depicted in Fig. [5j we use 
the Bogoliubov matrix (|98|) . Positive frequency "in" modes have amplitudes Ai„ = 1, Bi n = 0. The coefficients A out 
and B ou t are found by solving the system 



Aout i 1 ( ~l~ ^out ^out 1/1 



Bout J ly/QinQout V ^out ^in "f~ ^out J \ 



(99) 



whose solutions are 



* ~1~ ^out * j-, ^cmi ^ /inn\ 

Aout - , = = ttfefe , -Oout - ,„ „ = = -Pk-k • (1UUJ 

Z-y a ijj^a 'out ^"v ^ ^in^ ^out 
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FIG. 5: 'in' and 'out' basis in the temporal step-like discontinuity. 



These coefficients satisfy the unitarity condition 

\A out \ 2 - \B out \ 2 = K fc | 2 - |/3 fc _ fe | 2 = 1 , 



(101) 



where the minus sign means that the B ou t is associated to negative norm modes. 

Positive frequency "out" modes are characterized by A ou t — 1, B ou t — 0. The coefficients Ai n and Bi n are found 
by solving the system 



1 



fir. 



B, 



which gives 



fiin ~t~ fioui 



Bj. 



fi ; 



fi„ 



From these results, we see that the "in" and the "out" modes are related by the relations 

a ,in n * „.out q „,out* 

u k = a kk u k - Pk-kU_ k , 
and, considering the "in" and "out" decompositions of the field (/) 



in(out) 



dk 



in(out) in(out) 
k U k,<p 



a 



in(out)^ in(out)* 



k.^p 



we find the relation between the "in" and "out" set of operators, namely 



-out 
a k 



a* 

Pk-k a -k 



(102) 



(103) 



(104) 



(105) 



(106) 



The fact that both anhilitation and creation operators enter in the r.h.s. of the above equation means that the two 
decompositions (I105[) are inequivalent and that \in) ^ \out). 
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Density-density correlations 

The analysis of the density-density correlation is similar to the one performed in the hydrodynamic case, see [l]. We 
first write down the operator h 1 in the "out" decomposition 



h (t, x) — n 



dk 



and then use relation f|106D . For the two point function of n 1 in the |in) state we have 

/oo 
dk { [a% k (v%l + u$) - Pk-k(u X<f, + (*. *) x 

-OO 

x + - K_ k (u°_%t + u°_% v )] (f, x') + c.c.} 



(107) 



(108) 



This integral is well approximated by its hydrodynamical limit and the features of the density-density correlations 
are discussed in [26[ and [l|. 

B. Subsonic-supersonic configurations 

This case, which is relevant for the calculation of section [V] consists in a configuration made of an "in" subsonic 
region and an "out" supersonic one (ci n > \v\, c ou t < \v\). In the "in" region the analysis is the same as in the 
previous subsection. In the "out" (supersonic) one the dispersion relation (fT4|) shows new features with respect to 
the analysis in the hydrodynamic limit. From Fig. [S] we see that, for |fc| < |fc maa; |, the analysis is similar to that of 
the previous subsection, with the important difference that both modes are dragged by the flow and move to the left, 
whereas, when |fc| > |fc ma2; |, the supersonic modes &3(> 0) and k^{< 0) (in the language of subsection MI BP become 
able to propagate to the right upstream (from now on we find more convenient to work with positive fc, and indicate 
negative k with —k). The way in which the "in" modes propagate in the "out" region is shown in Fig. [6l These 
features become very important for the analysis of the temporal formation of acoustic black holes of section |Vl 

k <C k max 



-P-k 





k ^ kji 





FIG. 6: Evolution of 'in' modes for different values of k in the case of a supersonic 'out' region. 



For k > k max , an initial left- moving mode decomposes into a positive norm left- moving component plus a k^ 



negative norm one, with amplitudes A out and A™ respectively. These are found by solving 



1 



A ou t 

A4 J 2y/fli n Q out 



^out 
^out 



(109) 
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which yields the solutions 

A Qin + ^out * , ou t "in — "out _ Q /iim 

Aout = = a_ k _ k , A 4 = = . (110) 

These satisfy the unitarity condition 

\A out \ 2 - \A™*\ 2 = |a_ fe _ fc | 2 - |/3 fc _ fc | 2 = 1 . (Ill) 

An initial right-moving mode splits instead into a positive norm right moving k$ mode, with amplitude A^ ut plus a 
negative norm left moving one A ou t, which are found by solving 



A ut \ 1 j Qin + Qout ^in ~ ^out \ f 

A° 3 Ut J ~ 2y/n m il ou t \ ^in ~ Q out ^in + ^out J \ 1 



(112) 



the solutions are 



a out ^in ^out a A + ^out * ; 110 \ 

3 = /O. O = ~P~ kk ' ° ut = O /O- O = akk ■ ( 3) 

Eqs. (|110p and (111 3[) are the crucial formulas that we shall need in the next section to consider the temporal formation 
of acoustic black hole-like configurations. 

V. DENSITY-DENSITY CORRELATIONS IN THE FORMATION OF ACOUSTIC BLACK HOLE-LIKE 

CONFIGURATIONS 

In this section, with the help of the thorough analysis of the previous two sections, we will study the main Hawking 
signal in the more involved situation where an initial homogeneous subsonic flow turns supersonic in some region. We 
will model this situation with a temporal step-like discontinuity at t = (temporal formation) followed by a spatial 
step-like discontinuity at x — separating a subsonic and a supersonic region. The model we shall consider is sketched 
in Fig. where c r = Ci n . To study the propagation of modes solutions to Eqs. ([7]) for all x and t we need to impose 



t = 1 = 



FIG. 7: Temporal formation of a spatial step-like discontinuity temporally formed (c r = Ci„). 

matching conditions (IMt at t = at fixed k (only those for x < are non trivial), and then ((20)) at fixed to at x — 
(and t > 0). The behaviour at x — t = is more delicate because it depends the way we approach it. A detailed 
analysis of what happens for the case of subsonic flows was carried out in [l| by explicitly constructing the 'in' modes 
basis. As modes transiting through the origin only affect transient behaviours in the correlations patterns, in this 
section we will rather focus on those modes solutions which give the leading contribution to the main Hawking signal. 
We saw in the stationary analysis of section III.B that this is given by the evolution of the modes u^ m , u^ m * for 
w small and, consequently, > |A: mQ;! ;|. In turn, as shown in Fig. [5] such modes are generated by 'in' modes in the 
homogeneous t < region with the same value of k crossing the temporal step-like discontinuity on the x < side. 
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In our analysis we shall need to consider a transition from the k to the uj basis. The relations between modes and 
operators in the two basis are 



I duj 



-a k . 



(114) 



To construct the two point function (in|n 1 (t, x)n 1 (t' , x')|in) we proceed as usual by decomposing h 1 in the "out" uj 
basis 



h% x)= n ' ' k-tc 1 + v^T) + ar^ccr* + <r*) + ^wr* + <: ut ) + h - c -l - ( n5 ) 



and by relating the a™', a°"'^ operators to the a™, in the "in" (t < 0) region. This is done in two steps. First, the 
analysis in subsection IIII Bl provides for the relation between "out" and "in" uj basis in the t > region. In particular 
we have 



"-V.OUt 


= «? n - 























(116) 
(117) 
(118) 



From the values of the amplitudes in the above equation (given in subsection IIII Bh we see that the terms multiplying 
a^j m are subleading with respect to those multiplying d^ m and d^ m ^. 

Next, we need to jump from the "in" w-basis to the fc-basis needed to address the temporal step-like discontinuity. 
The relevant terms in h 1 in our analysis are 



n (t,x) 



dk-i 



(119) 



where = —k 3 . This is to be matched, at t = 0, at the relevant values of k, with the "in" decomposition (t < 0) 



n 1 ^, x) = 



dk 



f.inf in , in \ , -in (in , in \ , ~mf / m* , m* \ , ~mj (in* . in* 
a /c l u fc,0 + u /c, v J + a -k\ u -k,rp + u -k,ip) + a fc + + a -fc \ u -k,4, + u -fc, 



(120) 



The relations between the fc operators before and after the temporal step-like discontinuity are given by f| 106[) with 
k± = -k 3 : 



d k3 = a(k 3 )d k - P*(-k 3 )a?_ k , 
alfe 3 = -/3(-k 3 )a-k + a*{k 3 )a\ , 
where the Bogoliubov coefficients are given by (jllQI) and (|113[) 



a 



*2<\f Qin^out 



(121) 



(122) 



Here fl = \uj — vk\ is calculated before (f2j„) and after (f2 ou t) the temporal discontinuity. 

Let us now go back to the uj basis (the general relation between modes and operators in the uj and k basis is given 
in (11141) '). As shown in Fig. [SJ a fixed, positive value of uj corresponds to two values of k, namely £3 and k' A . We thus 
write 



n{t,x)= duj[a u (u U}<f> + u ) + a u '(^ +\ V )+^J- (123) 

Jo 

Defining k A — —k 3 and £4 = — k' 3 , the following properties are valid (we do not write explicitly the normalizations) 

u 3in* _ ^~iajt+ik 3 (ui)x^ _ ^iut—ik 3 (ui)x _ ^-i(—uj)t+i(-k 3 (uj))x _ e — i(-u})t+ik i (—tij)x _ u 4in* (124) 



and 



4m* — iujt+ik f 4 (ui)x — i{— u))t+i{ — fc 3 (— <^))x 3m* 



(125) 
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Therefore the density fluctuation operator n 1 turns into: 

/ j £ 3m / 3m , 3m \ , ^ 4m t / 4m* , 4-in* \ , * 4m t / 4m* , 4in* \ , 5 3m / 3m , 3m \ /-1 oc\ 

n(t,x)= / a w («^ + «^) + a, '(^ + u ^,v > + ( w -c^ + U-u,*) + «-a;("- w ^ + u ^^) -( 12 6) 



Since the w decomposition requires two values of k, the relation between the w operators and the k ones before the 
temporal discontinuity will involve relations (|121l) with different values of k, namely k 3 = k and k' A = —k': 



, / duj a a* l 1, ^ / dt ^at 



«( fc 3)\/ — afc 3 - P*(- k 3)y ^ u -k 3 ' 

„ « , , / , , , I duj 



dk 3 



a 4»t = _ /3( _ 4( _ a;)) J^i a _ fe , +a * (4( _ u;)) J^i a t / 



(127) 
(128) 



We compute now the Bogoliubov coefficients appearing above. Let us start with a(k 3 ) and /3( — /C3). By using again 
the fact that k is conserved in the temporal step-like discontinuity (fc = ^3) and the expression of £3 for small uj 

2\/v 2 — c 2 

(k 3 = v e;6 ' + -^j), which gives the main contribution to the density- density correlations, we can write fJ m and 



Qout = OJ — V 



2 _ r 1 



c/6 



n 



in — ^-in 



Cm A/ H - 



(129) 



(130) 



Notice that we cannot use the perturbative expressions in the "in" region, since here we are beyond the small frequency 
regime. Expanding up to uj we finally obtain: 



a(fc 3 ) 



-v + -cf + cfn , c /v /= ^(cf - cf n )(v + ^v 2 - cf + 



(131) 



2^{V 2 - C f+ Cfj 1 ^ 8v 2 ^TZ^f ( V 2 _ c 2 + c 2 n)5 /4 

s = V + V 7 ^ - c f + 4> , ClV^jcf ~ cj n )(-V + y^ 2 - C 2 + gjjj^ 

1 3 ' 2^( W 2_ c 2 +c 2jl/4 8w 2y^23^ (w 2_ c 2 + c 2j 5 /4 

Let us compute now a(k' 3 ) and /3(— k' 3 ). By using the fact that k is conserved in the temporal step-like discontinuity 
(— k = ki) and the expression of — k 3 for small uj (— fc 3 (— uj) = — — L + ^3^2) we have 



-uj + i; 



Q 2 



^1 = c in \ k 2 



u 2 — cf 



q6 



c 2 — v 2 



(132) 
(133) 
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By expanding (|122[) up to ui these expressions we finally obtain, at that perturbative level, 

a(k' 3 ) = a(k 3 ) , P(-k' 3 ) = . 



(134) 



We have now all the ingredients to calculate the main contribution to the Hawking signal in the density-density 
correlation for the temporally formed step. We study again the correlation between the modes u™ r and u™ 1 * . As at 
the end of subsection MI 51 x (x 1 ) is a point in the left (right) region. The two-point function reads 



m\{n\t,x),h\t' ,x')}\m){ul r ^ Ul 



ul* \ 



du) 



(<*<H^)| 2 + A*A L u *\p(-k 3 )\ 2 ) + <„)(*, x)(uj;; ut + u^frx')- 



+ (A^A^*\a(k 3 )\ 2 + A^*\P(-k' 3 )\ 2 ) (u u J : ; ut * + u^)(t )iB )«; + <£««*)(*, s') + 



The products of the amplitudes are related by 



r ur.out ul.out~\ 



o 



A L * A R 



A 1 * A r 



0, 



(135) 



(136) 



where we have used the relation between "in" and "out" ui operators. We neglect the subleading term A l *A r u and 
take into account that at leading order A l u *A r u is real. Thus, we find 



(M{hKt,x),h\if,x')}\m){ul r ^ul 



x A' u *A r u Re [« + <^)(i, + *')] H 

By taking in account that 

Hk 3 )\ 2 + \a{k' 3 )\ 2 + \p{-k 3 )\ 2 + \/3(-k' 3 )\ 2 



du { (\a(k 3 )\ 2 + \a(k' 3 )\ 2 + \P(-k 3 )\ 2 + m-k' 3 )\ 2 ) x 
c.}. 



2v 2 



v^/v 2 -cf+ 4 



(137) 



(138) 



we can finally write down the leading order contribution to G^ 2 \ namely 



G i2 \t-x,x') = 



{v 2 -c 2 fl 2 {c 2 



2v 2 



v v-\-c r 



v+ci ■ 



47m 2vci(v + ci)(v - c r )(c; - c r ) ^ v 2 - cf + c 2 



(139) 



v+ci 



which modifies the stationary correlation (|86p by the factor (|138l) that comes from the effect of the temporal formation. 
In Fig. (JSJ) we display the plots of Eq. (|139[) . and of the numerical counterpart along the direction x = x' — 1. The 
picture shows a good agreement, which confirms that the analytic approximation adopted in this paper is good enough 
to capture the essential features of the correlations. Good agreement exists also for different cuts, for completeness a 
3D contour plot is given in Fig. [TOl 

Finally, in Fig. ([lip we confront the signal between the eternal step and the temporally formed one. As one can 
see already from the analytic approximation, the temporal formation of the step yields an amplification of the signal. 



VI. FINAL COMMENTS 



In this paper we have studied in detail the formation of acoustic black hole-like configurations in BECs using step-like 
discontinuities. The Hawking signal in the stationary case ((56")) and in the case of temporal formation (| 139[) have 
stationary peaks (at = — f— ) of order 0(uj max ) ~ 0(l/£), which lie well inside the non perturbative regime in £. 

The results in the hydrodynamical limit of [18| showed instead a peak of order k 2 , where K = 4r\x=o is the surface 
gravity of the horizon. It is clear that in the approximation of spatial step-like discontinuities we are working with 
the surface gravity is formally infinite and therefore our expression (|139|) (and also (|86p ) regularizes the result of 18 1 
in the k — > oo limit, in agreement with the numerical results of 19]. 

In [l[, it was noted that a simple recipe to take into account a smooth transition region in c(x) around x = of 
width a x and surface gravity n ~ — is to introduce a cut-off of order n in the uj integral of (|137p by multiplying 

the integrand by the function e _w / K . The interplay between ui max and K is such that the final peak is of order 
k(1 — e^ u}max ^ K ), which has the correct k » Lo max limit, i.e. uj max . However, it is not able to make contact with the 
results of the hydrodymamic limit since, when uj max » n, we have a behavior in k which is linear and not quadratic. 
Thus, it would be interesting to find an analytical formula capable to interpolate successfully between these two limits. 
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FIG. 9: Comparison between the plots of Eq. (|139[) and of the numerical counterpart along the direction x — x' — 1, where x is in 
units of the length £. We adopted the following numerical values: v = —1.01, c; = — v/4, c r — —5v/3, n — 5.1, m = 20.1. 
With these choices £ ~ 0.03 with h = 1. These values have been chosen of the same order of the ones used in the simulations 
studied in 19|], so that they qualitatively match the results of [2C| |. 




FIG. 10: Comparison between the plot of Eq. (|139[) and of the numerical counterpart, obtained by numerically solving the 
integrals without truncating the expressions for the momenta as in appendix C. The variable y corresponds to x' . The values 
are as in Fig. ©. 



Acknowledgments 

We thank R. Balbinot, I. Carusotto, R. Parentani, and A. Recati for useful discussions. A. F. wishes to thank 
Generalitat Valenciana for financial support. A. F. and C. M. are supported in part by MICINN grant FIS2008- 
06078-C03-02. M. R. is supported by the Swiss National Science Foundation. 



25 




FIG. 11: Comparison between the (numerical) plots of the two-point function for the eternal step (|86[) and for the temporally 
formed step ()139[) along the line x = x' — 1 where x is in units of the length £ . The choice of the parameters is the same as in 
Fig. ©. 



Appendix A 

In this appendix we construct the exact "in" and "out" basis for the spatial step-like discontinuities at x = with 
v = (perturbative results in z\ = £w/c are given in the main text), which have a special interest for the validity of 
the unitarity relations. The scattering matrix determined by the junction conditions is given by 



Mscatt = W7 X Wr 



(Al) 



where Wi and W r are given by (|22l) and (1231) respectively. We recall that the structure of these matrices is determined 
uniquely by the matching conditions and by the solutions to the dispersion relation (|14p on the two sides. For v = 0, 
this equation reduces to 



2 2 

ur = c t 



(k l ) 



2 2 



r\2 , Cr(^ r ) 4 "' 



(«0 



(A2) 



where c ri ; are the speed of sound on the right-hand side and on the left-hand side of the step respectively. The 
solutions are, on the left-hand side 



(A3) 



K d,g 



± 



;V2 



» 2 % 



(A4) 



Similarly, on the right-hand side, we have 



K d,g 



= ± 



± 



V2 



-1 



1T\ 1 



2t2 



^ 2 e 



2t2 



UJ 2 t 



(A5) 



(A6) 
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As we have k v 
to 



-k,, and k d 



k„, and hence D u — D v (see Eq. (1151) '). Therefore, the matrices Wi and W r simplify 



Wi = 



and 



W r = 



( U u 


D l 

U u 




1*1 \ 




-*K D u 


" 7 / /"'/ 

~ lk d G <p 






h u 


fi 


r]l 

% 








lK g a V / 


( Dl 


D' u 


dl 


G } r \ 


iKD r u 




ikty - 


-ik g Gf 


K 


K 


d r v 






-iKK 


ik r d d r v ~ 


-ik r g G r v J 



(A7) 



(A8) 



Let us now construct explicitly the 'in' and 'out' modes (see details in section 3.1). 
• Mode 

Matching conditions at the step dictates that 



1 


V4/ 



= M m 




(A9) 



and R' = A l v and T = A r u are the reflection and transmission coefficients respectively. By solving the system we find 

(k v — fc£) 2 



\Rf 



4[(fc[) 2 -^) 2 ] [(K) 2 ~ (K) 2 ] \(K) 2 + 



2u 



I \2 



[(K) 2 - (k l d ) 2 ] [(K) 2 - (K) 2 ] ( k l) 2 + lt H 



I \2 



(Dl) 



(A10) 
(All) 



Note that both R' and T" do not depend on the normalizations d£ y To simplify the above expressions, we first use 
the definitions of D l u and D r u displayed in (IT5l) to find 



| T /|2 = 



^K) 2 [{k v ) 2 -{k d f] [{k l v f-{klf 



[(K) 2 ~ (k V d ) 2 ] [(K) 2 - (k 1 ^] [ki + kr] 

By calculating explicitly dk l v jdk r v and by using the identities 



dki 



dkl 



2 + c 2 



I \2 



y/uj 2 £ 2 + c? r = c r 



UK) 



r\2 



we find 



dk l v _ tfk r MrK) 2 + 2] 
dK ^l[(W+2] 



By noting further that 



we can write \T'\ as 



exit 



\T'\ 2 = 



±{Kf WI - (k r d ) 2 } 

K[(K) 2 - {k l d y} [K + K] 



(A12) 



(A13) 



(A14) 



(A15) 



(A16) 
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Finally, with the relations 



which can be easily proved with Eqs. (|A2[) and either (|A3I) or (IA5|) . we find that 



|Tf " «W ■ (A18 » 



Thus, \T'\ 2 + \R'\ 2 = 1 



• Mode u v J% 



The scattering matrix is still given by (|A1|) , but now the system to solve is 



( Al A 




\ Al J 




(A19) 



with R = A r u and T — A l v . Despite the fact that the system is different, we find the same results as in the previous 
case, namely 



\R\ 2 = 



\T\ 



4[(fc;) 2 - (fc-) 2 ] [(fc^) 2 - [k l d f] [{k[f + f^] 2 (*g 

[(fc[) 2 - m 2 ] my - (k l d y] \{ky + ^1 2 [fci + / 



(D r J 



(Dl) 



(A20) 
(A21) 



thus even in this case unitarity holds. It is easy to see that these expressions can be obtained by the ones in the case 
u u,m ky swa pp m g r -f> i, so that the proof that |i?| 2 + |T| 2 = 1 follows immediately. The construction of w"'^ 1 and 



u "j°(p t 1S now straightforward. 



28 



Mode <>7 



For u"'T* we solve 

u),<p 



( °A 



\ Al J 



where we identify R* = A\ and T* = A\ and find 



= M S , 



V o 



liTI 2 

| T .|2 



4[«) 2 -(fcS) 2 ][«) 2 -(^) 2 ] 



(fc,) 2 + FT- 



m) 2 -{k r d ) 2 Wv) 2 -W] (K) 2 + ft [K + K 



• Mode u 



UJ,(f> 



For u v J we solve 



( ^ 

V 





s 



( K 
o 

, AS 
V o 



where R'* = A l u and T'* = A r v and get 

(k l v — k T v ) 2 



,/* 1 2 



| T /*|2 



4 t(fc[) 2 - [(Kr - (fc S ) 2 ] \(k? + ^t]\k{f 



[(K) 2 - (k l d ) 2 ] [{KY - (fcS) 2 ] [«) 2 + {%] [*i + fcl] 

In both cases unitarity relations are satisfied. 



Ml! 
(i?r) 2 



(A22) 



(A23) 
(A24) 



(A25) 



(A26) 
(A27) 
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Appendix B 

In this appendix we give the perturbative results for the construction of the "in" modes for spatial step-like discon- 
tinuities at x — for v =/= ( in the subsonic-subsonic case). The details of how to construct them are given in 
subsection 1111 Al with the help of Fig. 1. For simplicity, we will only give explicitly the amplitudes of the propagating 
it, v modes. 

i 

l U,<f> 

By solving the system (|25|) we find, for the propagating modes, at O(zf) (where zi = £zw/cz) 
, _ 2^/cWr . 4' 2 ( c l - <V) c r (\/ c t -v 2 - V c l - v2 ) z l 

Cl+Cr ( C J -V)(V- C r ) (cj + C r ) yj ' (v 2 - cf) C r (v 2 - C 2 ) 

Zfc^id-Cr) 2 r _ v 3 {v _ Cr) 3 iv + Cr) 2 + 



8 (V ~ Cl) 3 (V + C,) 2 (V - C r ) *Cr /2 (V + C r ) 2 (q + C r ) 
+ V 2 Ci (V - C r ) 3 + C r ) 2 + cf (v 3 - V 2 C r + 3vcf + C 3 ,) + Vcf (~V 3 + C r (v 2 + C r {v + 3c r ))) + 

+ cf (-2v 5 + c r [2v 4 + c r {v + c r ) ( -5v 2 + 2vc r + c 2 r - aJ(v 2 -c 2 )(v 2 -c 2 )) )) + 



+ vcf ^2v 5 + c r ^-2w 4 + c r (v + c r ) (~3v 2 - 2vc r + 3c 2 . - iyj(v 2 - cf) (v 2 - C 2 )^j ^ } 



(Bl) 



Cf (Cl - C r ) C r (v 2 - C 2 + y/ (cf - V 2 ) ( C 2 - V 2 )) 



+ 



Cl + C r yjcf - V 2 (Cl + C r ) (v 2 - C 2 ) 2 

cf (q - c r ) z 2 



+ 



r, — — 9 ^ [2cf C 3 - W 2 (« 2 " C 2 ) 3 + Cf (-U 4 + cf ) + 

4(w 2 -c 2 ) 2 c I .(Q+c r )(i; 2 -c 2 ) 3 L v ; v ; 

+ c 2 (2« 6 - v A c 2 r - 2v 2 c A r + cl) + 2c 3 c 3 (-3« 2 + c 2 - 2yJ(v 2 - cf) (v 2 - c 2 )) + 



+ 2v 2 cicl ( -c 2 + 2[v 2 + J(v 2 - cf) (v 2 - cl) 



(B2) 



The important check is the unitarity relation \A V \ + = 1, which is satisfied quite non-trivially at O(zf), as 

l |2 _ 4c;C r zfcf (Cl - C r ) 2 



1-4 



\K 



(Cl + C r ) 2 2(v- Cl) 3 (V + Cl) 2 (V - C r ) 3 C r (v + C r ) 2 (cj + C r ) 
U 2 Q (« - Cr) 3 + C r ) 2 + vcf (v - C r ) (v 2 + cf.) - cf (v - C r ) (v 2 + cf) + 

cf(ci - V) (2v 5 + C r (-2v i + C r (v + C r ) (v - C r ) 2 )) , 

Ci - Cr ° 2 ( Q _ Cr ) ° r ~°r + "S^ 2 ~ ^ 2 _ C r)) ^ 

Cl + C r y/v 2 - cf (Cj + C r ) {v 2 - C 2 ) 2 

cf (Cj -c r )z 2 



2 [w 3 (t! - C r ) 3 (V + C r ) 2 + 



(B3) 



4 (?; 2 - cf ) 2 c r (cj + c r ) (w 2 - c 2 ) 



^[2cfc 3 ~, 2 (, 2 -c 2 ) 3 + cfH 4 + cf) + 



(2v 6 - « 4 c 2 - 2t; 2 cf + 4) + 2c 3 c 3 ^3« 2 + c 2 - 2y / (« 2 - c 2 ) ( V 2 - c 2 )) + 



+ 2 V 2 C/ cf -c 2 + 2 U 2 + ./(w 2 - c 2 ) (w 2 - c 2 ) 



(B4) 
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• Modes it™" 



The construction proceeds from Eq (f3"2"j) , which gives, at O(zf) 



Or - Q . 1 
I 



(ci - Or) [J c'f - V 2 - ^Jc 2 r ~ V 2 ^ Z\ 



Cl+C r ( c 2 _ „2) 3/2 ( Q + ^5T 



cf (Q - Or) zf 



V (v 2 -c 2 r ) 2 + v A cf (3v 2 - c 2 r ) + 2c\c r {c 2 , - v 2 ) + 



A., 



4 (cf - v 2 ) 3 c r (ci + Or) (v 2 - c 2 ) 2 
+ cf {v 2 + c 2 r ) + cf (-3 V 4 - 2v 2 c 2 r + cf) + 2c 3 c r (c 2 r - v 2 ) [c 2 r - 2 (V + yf(v* - cf) (v 2 - c 2 )) ) 

(ci - Or) {^Jcf - V 2 - ^C 2 - V 2 ^j Z\ 



2y/ClC r 

Cl+ Or (v + Ci) 3 / 2 (v + Or) (q + C r ) y/c t ( C/ - v) (c 2 - V 2 ) 



+ 



cf /2 (q- c r ) 2 zf 



[v 3 (V ~ C r ) 2 (V + C r ) 3 + 



(v - Ci) 2 (v + Ci) 3 (v - C r ) 2 cV 2 [v + Or) 3 (c ; + C r ) 

+ v 2 c t (v - Or) 2 (v + Or) 3 + cf (v 3 + c r (v 2 + 3vc r - c 2 )) + ucf (u 3 + c r (v 2 - vc r + 3c 2 .)) + 



- vcl \ 2v b + c r i^v* -(v- Or) Or {3v 2 - 2vc r - ic 2 r + 4^/ (v 2 - cj) {v 2 - c 2 )jj ) 

- cf (2v 5 + c T \2v 4 + (v - c r ) c r (^5v 2 + 2vc r -c 2 r + 4yV - cf) {v 2 - c 2 )^ J J 
The unitarity relation \A l v \ 2 + |A£| 2 = 1 is again non trivially satisfied at O(zf), being 



C T - Cl 



2 Z 2 cf (Cj - Or) 2 



A I 2 2\ 2 1 42/2 o2\, 

" (v — c r ) + v Ci [c r — 3v ) + 



ci+c r J 2 (v 2 - cf) 3 c r (ci + Or) 2 (v 2 - cl) 2 
cf (v 2 +c 2 r )+ct(3v 4 + 2v 2 cl-c r )], 

AciCr zfcf (Cl - Or) 2 



[v 3 (v-C r ) 2 (v + Cr) 3 + 



(C; + C r ) 2 2 (v - Cl) 2 (v + Cl) 3 (v - C r ) 2 C r (v + C r ) 3 (q + C r ) ' 
V 2 Ci (V - C r ) 2 (V + C r ) 3 + Vcf {v + C r ) (v 2 + C 2 ) + cf (v + C r ) (v 2 + C 2 ) + 

cf(ci + V) (-2V 5 + C r (^-2v A - (ll - C r ) C r (v + C r ) 2 )) 



(B5) 



(B6) 



(B7) 



(B8) 



Appendix C 

In this appendix we extend the leading order results of the calculations of the amplitudes of the propagating modes 
w w™ an d m the case of the subsonic-supersonic spatial step-like discontinuity. With these, we are able to check 

the unitarity relations. 

.4m* 

The amplitudes depicted in Fig. 4 are: 

/2c7(w 2 -c 2 ) 3 /> + c r ) 



A - - 7^M^{^-'^-v +(A+w) ^' (C1) 

K = J V lZ'* f ' >+ r ^ = (-M^ - ij^cf) + (C + iD)^i , (C2) 

C{ y/2zi(ci + C r )y/cl- V z \ J 

Ai - - ,v ^) + {E + iF) ^' ■ (C3) 
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where 



.4 



(-8w 6 + 2vcf (-2v + cv) + 2w 2 c r (v 3 + 4v 2 c r - 2c 3 ,) + cf (7v 4 ~ c r (Av 3 - 2v 2 c r + c 3 ))) c 4 ^c7zi 



2V2v (v 2 - cf) 3 / 4 (v - Cr) 2 (« + cv) (c 2 - C 2) 



B = - 
C 



(8v 3 + 6v 2 c r + cf (-V + Cr)) Ciy/c r (-V 2 + C 2 ) Zi 



2V2V (V 2 - C 2 ) V4 („2 _ c 2) ( c 2 _ c 2) 

1 



8u" + 2v 2 c\ + 2vcic 2 r (v 2 - c z r ) + 2vct(-v 2 + c z r ) + 



2^„.2 2\ , r,„.„3/ 2 , 2\ 



L> = 

E = 



AV2v(v 2 - C 2 ) 3 /4( v - C r ) 2 (v + Cr){ci + C r ) 

+ cf(-7v 4 + 2v 3 c r + 2v 2 c 2 r - 2vc 3 r + cf) + 2v 2 c r (-v 3 + c r (-Av 2 + c r (v + c r ))) 

(8« 3 + 6v 2 C r + cf (— U + C r )) y/ClZl 



'ClZl 



iV2v {V 2 - Cf) V4 (q + Cr ) ^(- V 2 +c 2) 

1 



8v 6 + v 2 cf (~7v 2 + 2c t (v + c/)) + 



F 



AV2V (V 2 - cf) 3 / 4 (V - C r ) 2 ( C j - Cr) (v + C r ) 

+ 2v 3 (-v 2 + cf) c r - 2v (Av 3 + a (v 2 - vci + cf)) c 2 + 2v (v 2 - cf) c 3 + (2v 2 + 2vq + cf) 4 
(8v 3 + 6v 2 c r + cf (—v + c r )) y/cjzl 



'ClZl 



AV2v (v 2 - cf) V4 ( Q _ Cr ) V /(Z^2 + c 2) ■ 
Their squared modulus read 

2c r {v 2 -cf) 3 ' 2 {v + c r ) 



1^2 I 2 



G 



V |2 



cfzi(cf - C 2 )(v - C r ) 
[v 2_ c 2 ) V2 {v + Cr){Cr _ Ci) 

2cfzi(c l + c r )(c r - v) 



H 



where 



G 
H 



|Mi)n<oi2 = {v 2 - cf) 3/2 {v + c r ){g + c r ) | ] . 

2 2cfzi(ci - C r )(v - C r ) 



(-2v i + Av 2 c 2 r + cf (-Sv 2 + c 2 r )) 



V (v - C r ) 2 (-cf + C 2 ) 

(c t - Cr) (2v 4 - 2v 2 c r + cv) + cf (v 2 + 2vc r - cf) + 2vq (-v 2 + cf)) 

AVCI (v - C r ) 2 (cj + Cr) 

I _ (c t + c r ) (2v 4 - 2v 2 c r (v + c r ) + 2vq (v 2 - cf) + cf (v 2 + 2vc r - cf)) 

AvCl (v - C r ) 2 (ci - C r ) 

The above amplitudes satisfy the unitarity condition \A r u \ 2 + \A l v \ 2 — \A l u \ 2 = — 1 at this perturbative level. 



(C4) 

(C5) 
(C6) 

(C7) 
(C8) 

(C9) 

(CIO) 
(Cll) 
(C12) 

(C13) 
(C14) 
(C15) 



• Modes u 



3,m 



The amplitudes sketched in Fig. [4] turn out to be 

y/2c~ r (v 2 -Cf) 3 ' 4 {v + C r ) 
Cly/zl(cl - cf)y/c 2 - V 2 
(V 2 -Cf) 3 / 4 (V + C r ) 



A' 



\Jc 2 r - V 2 + iJv 2 -cf) - {A- iB)y/zi , 



C^ 2 yj2zi(ci + C r )^/c 2 r - V 2 

(v 2 -cf) 3 l 4 {v + c r ) 

C^ 2 y/2zi(ci - Cr)y/(% - V 2 



^cj^ + iJv 2 - cf] -{E- iF)^zl , 



(C16) 
(C17) 
(C18) 



.32 



where 



i^i 2 



2c r (v 2 -c 2 f/ 2 (v + c r ) c 

cfzi{cf - c 2 )(v - c r ) 
(v 2 -^fl 2 (v + c r )(c r - Cl ) 



H , 



(C19) 



(C20) 



\A L U \ 



2cfzi(ci + c. r )(c r - v) 
(v 2 -cf)^ 2 (v + c r )( Cl +c r ) 



I . 



(C21) 



2cfzj(Cj - Cr)(v - Cr) 



Again, one sees that \A^\ 2 + \A%\ 2 - \A^\ 2 = 1 is satisfied. 
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